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Abstract 

We introduce a model of self-repelling random walks where the 
short-range interaction between two elements of the chain decreases 
as a power of the difference in proper time. Analytic results on the 
exponent v are obtained. They are in good agreement with Monte 
Carlo simulations in two dimensions. A numerical study of the scaling 
functions and of the efficiency of the algorithm is also presented. 



1 Introduction 



Self-avoiding walks are a fascinating subject. Although they are much simpler 
than an Ising or Heisenberg spin system, they behave quite similarly in the 
critical region where the coherence length is much larger than the lattice 
spacing. There are deep reasons for this similarity. Indeed, using Symanzik 
representation M, |2[], the free energy of a spin system may be written as sum 
over partition functions of self-avoiding walks. Moreover the self-avoiding 
walks are exactly the limit of an 0(n) -invariant spin system for n going to 
zero 1,1,1,1. 

A long time ago Flory |7|, using simple minded approximations, found 
that the critical exponent v in dimension D less than or equal to 4, was 
given by the simple expression 

i/ = 2/(3 + £>) (1) 

This result is quite puzzling; it is exact for D =1, 2 and 4, but it is definitely 
wrong in D = 3 0, [10], [II]]. Moreover in 4 — e dimensions the exact results 
and Flory's one differ at first order in e f3|. The Flory approximation was 
strongly criticized by des Cloiseaux [ IS , |K| , but he found the puzzling result 



that a more accurate variational approximation leads in three dimensions to 
the rather bad value of v = 2/3. 

Quite recently it was found that des Cloiseaux' method can be extended 



to the case of polymers with long-range repulsion [ITJj]. Also in this case 
Flory's and des Cloiseaux' approaches give different results, however here 
it is possible to give very strong arguments suggesting that the variational 
approach of des Cloiseaux gives the exact result in the case where the space 
dimension becomes infinite. 

The putative exact infinite-dimensional results can be extrapolated at 
finite dimensions using arguments based on the renormalization group and 
performing the approximation directly at finite dimension. One finally finds 
that for an interaction which decreases with the distance as a power law with 



exponent A, v is given by the generalization of des Cloiseaux' formula |T4 | 

for A < 2 

for 2 < A < 4 (2) 
for A > 4 

with logarithmic corrections when A = 2, 4. 
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Some of these predictions have been numerically tested ||15|| . 

A very different case concerns a potential which is short-range in space but 
in which the interaction among different elements of the polymer decreases as 
a power of the distance along the chain. An interesting theoretical analysis 
of this model can be done and the predictions can be tested using numerical 
methods due to the short-range nature of the Hamiltonian. 

After this introduction we present the model in section 2. An approximate 
computation of the exponent v is presented in section 3. The algorithm we 
use for numerical simulations is described in section 4, while the results we 
obtain are presented in section 5. Some conclusions are presented in section 
6. 



2 The model 

Let C be the ensemble of random walks on the hypercubic lattice TP of N 
steps, starting from the origin. Thus, if uj G C 

u = {u ,ui,u 2 , ■ ■ ■ ,3n} (3) 

with uj q = the origin of the lattice, u)i G IP the location at time i of the 
walk and \uji — uJi-i\ = 1 for i — 1, . . . , N. 
We consider an interaction of the form 

1 N S- - 

Z i,j=0 \ l J\ 

This means that in the ensemble average (•) in C we associate to each walk 
the statistical weight 

e --ff/[w] 

m \u} = y r 1 (5) 

Zn[9\ 

where 

^] = E"" ff,H (6) 

ujGC 

To study the conformation of the walks in this ensemble we shall consider 
the square end-to-end distance 

RIM = (7) 
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and the square radius of gyration 




Both of them are believed to have the asymptotic behaviour 

(R 2 ) ~ N 2v (9) 

as N — > oo, with the same critical exponent v. 
We shall also consider the universal ratio 

For the ORW we have A = 1/6 ~ 0.16667, while for the SAW A is known 
only numerically and depends on the dimension of the embedding space. 
In two dimensions the most precise estimate is obtained by a Monte Carlo 
simulation |16[] whose result is A = 0.14026 ± 0.00011, where the error bar is 
95% level of confidence. 



3 Heuristic analysis 



In this Section we want to derive an estimate of the critical exponent v with 
a variational approach. Let us firstly discuss the case 6 



M 



Following Des Cloiseaux |TI| one considers random rings instead of ran- 
dom walks, as this fact does not change the value of v and simplifies the 
calculations by permitting Fourier analysis along the chain. In addition it 
is simpler to work in continuum space rather than on a lattice. In this case 
one must somehow regularize the 5-function appearing in (^). We will thus 
consider as equilibrium probability measure for the model 



dm(u) 



exp(— if ) d D uJid D uj 2 ■ ■ ■ d D uj^ 



(11) 



where the Hamiltonian is given by 



N 



2 H< 

1=0 



1 N 
,u)i - Qi-if + - g 

Z i,j=0 



V 


(UJi -ujj) 2 




i-j 


X 



(12) 
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Here V(x 2 ) can be any arbitrary short-range potential and for definiteness 
we will assume 

^ = exp (~h) 



V(x 2 



(13) 



The mean field approximation is based on a variational approach with a 

exp(— Hq) d D uJid D uj2 ■ ■ ■ d D UN 



Gaussian trial measure [12[ [14| 

1 



dm (u) 



with 



and 



N 



Hn 



Z ij=0 



(2 



7T 



,(iV+l)D/2 



(det G) 



D/2 



The function Gy is determined by minimizing the functional 

F[G] = (H-H ) -log Z 



(14) 

(15) 
(16) 
(17) 



where ( ) denotes the expectation value with respect to dm . 

Because of invariance under translations along the chain we have Gy = 
G(i — j) and by symmetry G(n) = G(—n). It is also convenient to introduce 
the Fourier transform 



N-l 



G(P) = E G{r)e 



IpT 



(18) 



T=0 



where p — 0, 2n/N, . . . , 27r(iV — 1)/N, and G(p) is real and positive due to 
the positive definiteness of G^-. 

The functional F can be computed as 



F[G] 



" ^G(p) (1 - cosp) - £ [1 + log(27r)] - £ f ^ logG(p) 

-7r Z7T Z Z J-7T Z7T 



i + 2 r ^ G(p)(1 _ cospr) 

a J-7T Z7T 



_> 



(19) 



where we have replaced sums over p with the corresponding integrals as we 
are interested in the regime of very large N. 
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The minimization condition becomes 
1 



G(p) 



= 2 (1 — cosp) 
g 00 1 

~X)~r (l-cospr) 



(20) 



T = l 



T' 



a j-tt 27r 



-, -D+2 
2 



The exponent v can be extracted, as easily seen, from the low-momentum 
behaviour of G{p) as 

G(P) ~ ^4tt (21) 



The analysis of ( p0[) is quite subtle and can be done following the original 
paper by Des Cloiseaux [12] and predicts v = 1/2 for D > 4 while for 
2 < D < A gives the following result |T4| 



v MF (\,D) 



1/2 

1 f2 - A) for A < ±(4 



for A > i(4-L>) 



For A = 0, D = 2 additional logarithmic corrections appear as 

(R 2 g ) ~N 2 /\ogN 
and analogously for A = |(4 — D) we get 

(^)~iV(logiV) 2 / B 



(22) 



(23) 



(24) 



Thus we obtain an ORW for A large enough. This has indeed to be expected, 
as we have the trivial rigorous bound 



H^u) < CN 



2-A 



(25) 



showing that at least for A > 2 the model is an ORW. 

On the other hand the mean-field result Vmf is definitely wrong for A = 0, 
D < 4, which corresponds to the SAW limit Jl7|, [ll|. In two and three di- 
mensions v S aw is equal respectively to 3/4 and approximately 3/5 (see ([!])) 
while vmf{Q) corresponds to 1 and 2/3. For D = 4 the mean-field approach 
correctly predicts logarithmic corrections to the random-walk behaviour, but 
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the power of the logarithm is higher than that obtained using the renormal- 
ization group which predicts |18| 



(i§~iV(logiV) 1/4 (26) 

Of course the variational approach will overestimate v also for small A. We do 
not have any theoretical control of this regime. However as we will discuss 
in Section 5, our numerical results are in reasonable agreement with the 
following conjecture for u(X, D) 

u(X, D) = min (u SAW , u MF (X, D)) (27) 



Let us finally notice that for large N, the mean-field equation (|20| ) has a 
general scaling solution. Indeed in the limit of small p, which corresponds to 
large N we can rewrite (|20|) as 



' P 2 (28) 



G(p) 



-g ^ P 2 r drr^ - r ^G( P ') (i - cospv 

2a Ji [a J-tt 2n 

It is easy to check that the general solution has the form 

G( P ) = -Jx(9P X - 2+D/2 ) (29) 
p A 

which gives 

{R 2 g )=Nf x (gN 2 - x - D l 2 ) (30) 

Let us now discuss the case 5^0. For simplicity we shall set since now on 
D = 2, the generalization to arbitrary D being trivial. 

Let us firstly consider the case in which 5 < 0, namely the case in which 
we are weakening the coupling constant when approaching the asymptotic 
limit. One can think at the model in the intermediate region in a perturbative 
expansion around the ORW, that is around g = 0. To evaluate the dimension 
of the coupling constant, remark that N 2 ~ uD is the asymptotic behaviour for 
large iV of the average number of intersections of two walks of iV steps and 
Hausdorff dimension \ jv on a lattice of dimension D. Then 

[ g ] = -[N](8- X + 2-uqrwD) = 2 + 25 -2X (31) 



D+2 
2 
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and 

[N] = — = -2 (32) 

v ORW 

Thus in the scaling region we expect that for gN s small 

(R 2 g ) = NF x (g,N) = Nf x (gN 1+5 - x ) (33) 

Let us notice that for 5 = this expression coincides with the scaling formula 
(j30|) . For g = the model is an ORW and thus we get f x (0) ^ 0. It follows 
that, when 5 < A — 1 (R 2 ) scales as N . On the other hand when A — 1 < 5 < 
the argument of the scaling function fx goes to infinity. Assuming in this 
limit f\(x) ~ x 13 we get 

(R 2 g ) ~ N l+ ^ 1+5 - x ^ (34) 

The value of (3 is computed using the conjectured value for v when 5 = 0. In 
this way we obtain that f\(x) scales for large x as 

x i-> ~x 2 d- A ) when < A < 1/2 /ocN 

MA)-i _ , ~ 7 , ( 35 ) 

x !- A ~x whenl/2<A<l 
Then we obtain the prediction 

f at§+2(T^a) when < A < 1/2, 5 < A - 1 
(^s) ~ | AT2+5-A when 1/2 < A < 1, 5<\-\ (36) 
I AT when A > 1 and A < 1, 5 > A - 1 

In the opposite situation (5 > 0) the coupling constant goes to infinity while 
approaching the asymptotic limit and we can expect that in this case the 
scaling behaviour is controlled by the SAW fixed point. Then 

[N] = = ~\ (37) 

and 

\g] = -[N]{5-\ + 2-u S A W D) (38) 
so that we expect that for gN s large 

(R 2 g ) = N 2 ^ F x ( 9l N) = Nlf x (gN^ +s - x ) (39) 
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Since for g going to infinity at fixed N we expect the model to describe a 
SAW, fx(x) must converge to a constant for large x. It immediately follows 
that for 5 > A — \ the behaviour is that of a SAW. On the other hand, when 
5 < A — 1/2 and A > 1/2, the argument of the scaling function goes to zero 
when iV goes to infinity. As for g = we have an ORW, we must have 
/a(0) = 0. Then, assuming for small x the scaling form f\(x) ~ x 13 we get 

(R 2 g ) ~ iV3/2+/3(l/2 + 5-A) (4Q) 

This exponent (3 is computed by requiring that for 5 = this formula repro- 
duces the result (|27|) . We thus get that the scaling function fx must scale for 
small argument as 

{ 4,{ v (\)~v SAW ) 
x i-» ~x whenl/2<A<l m 

x !- 2A ~ x 2 ^- 1 when A > 1 

Then for positive 5 we get 

f iV 3 / 2 when < A < 1/2 and A > 1/2, 5 > A - 1/2 
{R 2 g ) ~ I N 2+5 - x when 1/2 < A < 1, 5 < A - 1/2 (42) 
[jV 1+2 ^r when A > 1 5 < A - 1/2 

Let us notice that all these scaling arguments do not take into account loga- 
rithmic corrections which we expect to be present for those values of 5 where 
there is the transition to the purely SAW or ORW behaviour. 



4 The algorithm 



We have simulated the model by using the so-called pivot algorithm 19, 20 



2IJ , which is known to be extremely efficient for the simulation of SAWs with 
fixed number of steps and free end-points as the computer time necessary to 
produce an independent walk is of the order of the number of steps in the 
walk, which is also the best possible behaviour because this is the order of 
time necessary simply to write down all the steps. 

The algorithm is defined as follows [2T]. Choose at random a point along 



the walk (the pivot), but not the first or the last one. Let it be the k-th point 
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Uk, with < k < N. Then choose at random an element g in the symmetry 
group of the lattice and propose a new walk uj' defined by 

^ f uji for < i < k . . 

~ 1 o3 fc + - dJ fc ) for jfe + 1 < % < N 1 ' 

The new walk is accepted according to a Metropolis test in order to generate 
the desired statistical ensemble. 



For the limiting case of ORW and SAW it is known |21| that for a global 
observable A the integrated autocorrelation time r int) A scales for large number 
of steps according to 

Tint,A ~ N" (44) 

where (the exponent p should be the same for all global variables) 

v = 1 for 0RW (AS\ 

y 1 0.194 ±0.002 for SAW 1 ; 

For our models, as it can be seen from Table [1], where we report the integrated 
autocorrelation time for the end-to-end distance, we obtain similar results. 
When the values of the parameters are such that the walks are in the same 
universality class of ORW or SAW, the dynamical behaviour is compatible 
with (f45D- In the intermediate cases for which 1/2 < v < 3/4 the dynamic 
critical exponent p ranges correspondingly within the interval [0,0.19], and 
we find that it is in reasonable agreement with a linear interpolation of the 
form 

p v = 0.39 (2 v - 1) (46) 

Let us now come to the computational complexity. In the practical imple- 
mentation we used a hash-table with linear probing in order to check for 



the self-intersections of the walk. In the case of SAWs Madras and Sokal pi 
showed that it is particularly convenient to insert the points in the hash-table 
starting from the pivot point and working outward (and of course stopping 
as soon as an intersection is detected). In this way the mean work per move 
turns out to be of order iV 1-p , thus smaller than the work done by inserting 
the points without a special order which is of order N. We used in our case 
a similar trick. Indeed, chosen a random number r uniformly distributed 
in the unit interval, according to the Metropolis prescription, the proposed 
walk u' is accepted if 

H[lu']< S = H[uj]-\nr (47) 
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Our implementation works as follows: after the choice of the pivot point and 
of the transformation g we choose the random number r and compute the 
quantity S. Then we begin to construct the new walk and we insert the 
points in the hash-table starting from the pivot point and working outward. 
Whenever we find an intersection we sum to the energy of the new walk the 
contribution from that intersection, but we stop if the accumulated value 
is already larger than S, because the proposed walk has to be rejected Q. 
In the limiting case of very strong repulsion, when the universality class is 
the same of SAW, this implementation works exactly as the original one 
devised by Madras and Sokal. In the opposite limiting case, the ORW-one, 
practically all proposed walks will be accepted and thus no improvement can 
be expected. In the general situation we expect (although we didn't check) 
that the computational work scales as N^ 1 "^ where p is the dynamic critical 
exponent for global observables and thus, in all cases, that the computer time 
necessary to produce a statistically independent measurement is of order N. 



5 Numerical results 

In order to test the ideas presented in Sect. 3 we have performed an extensive 
Monte Carlo simulation on a square lattice on walks of lengths ranging from 
100 to 8000. The total CPU time for these runs was roughly 2000 hours of a 
VAX 6000-520. 

In the numerical simulation the first problem one has to deal with is the 
initialization, that is, for each value of A, 5, g and N one has to generate a 
starting walk. This was done in two different ways according to the value 
of N. When N was less than or equal to 2000 we generated a SAW using a 



dimerization routine [22, 23, El En. When instead N = 4000, 8000, as the 



dimerization routine is too costly (the computer time needed to generate a 
walk scales as r ~ N alog2 + , with a ~ 0.17 and b « 0.72 in two dimensions), 
we used the scanning method \25\ with scanning parameter equal to 3 None 

1 This procedure works because the interaction is repulsive. For attractive interactions, 
as for instance for SAWs or trails at the 0-point, the whole new walk must be defined and 
therefore the mean work per pivot-move will scale as the number of step N. 

2 The only reason why we used the dimerization routine for the low-iV runs was its 
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A 


0.00 


0.00 


0.00 


0.00 


0.25 


5 


-1.00 


-0.75 


-0.50 


-0.25 


0.00 


9 


1.00 


2.00 


3.00 


4.00 


1.00 


A 


= 100 


6.27(5) 


6.22(5) 


7.00(6) 


8.07(8) 


7.48(7) 


A 


= 200 


6.31(8) 


6.24(7) 


7.76(10) 


9.07(13) 


8.53(12) 


A 


= 500 


6.54(8) 


6.44(7) 


7.99(11) 


10.77(17) 


9.87(15) 


A 


= 1000 


6.53(8) 


6.61(8) 


8.62(12) 


11.63(20) 


11.19(18) 


A 


= 2000 


6.51(8) 


6.54(8) 


8.74(12) 


13.24(23) 


13.14(23) 


A 


= 4000 


6.40(8) 


6.64(9) 


9.59(14) 


14.71(27) 


15.37(29) 


Pu 


0.00 


0.05 


0.09 


0.14 


0.19 


p 


0.00(1) 


0.01(1) 


0.08(1) 


0.15(1) 


0.20(1) 



Table 1: Integrated autocorrelation times for the end-to-end distance for 
selected values of the parameters A, 5 and g and of the number of steps A 
of the walks. p u is the value obtained through a linear interpolation for the 
dynamic critical exponent, and p is the estimate based on a fit from our data. 



of these two methods generates a random sample of walks with the correct 
equilibrium distribution, although the first is exact in the SAW-limit and the 
second in the ORW-one. It was thus necessary to run a certain number of 
thermalization iterations before measuring. As the convergence to equilib- 
rium of a Markov chain is controlled by the exponential autocorrelation time 
r exp it is necessary to run a few r exp iterations to reach equilibrium. For the 
pivot algorithm r exp is proportional to A. For this reason we ran approxi- 
mately 10 AT pivot iterations for thermalization before measuring. For each 
value of A, 5, g and A we have then performed 10 6 iterations, except when 
A = 100 in which case the runs consist of 2 10 6 iterations. The integrated 
autocorrelation times for the squared end-to-end distance and for the square 
radius of gyration range from 3 to 20 and from 5 to 60 respectively. We have 
reported a few of them in Table 1. 

In Table [I] and Table |] we report the results of our runs for (R 2 e ) for dif- 
ferent values of the parameters, respectively for negative and positive values 

availability at the time of the runs. In retrospect it would have probably been better to 
use the scanning program in all cases. 
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of 6. 

We have performed least-squares regressions on these data in order to 
determine the critical exponent v. We fit (R 2 ,) to the Ansatz aN 2v by per- 
forming a weighted least-squares regression of its logarithm against logiV, 
using the a priori error bars on the raw data points to determine both the 
weights and the error bars. In order to control the systematic error due to 
corrections to scaling we have done various fits in which the data points with 
lowest iV were discarded. The results of these fits are reported in Table 4 
and Table 5 for various N cut , where N cut is the minimum N included in the 
fit. In many cases it is evident a systematic drift of the estimated exponent 
with N cut , an indication of strong corrections to the scaling. This effect is 
however strongly dependent on the value of g. When the expected value of 
v is different from 1/2 and 3/4 one observes that for g small the estimated 
value of v increases with N cut , while for g large the estimate decreases. For 
the intermediate values of g there is a flatness region where v remains approx- 
imately constant meaning that for these values the (/-dependent corrections 
are small compared to our statistical error. For every A and S we have made 
various runs for different values of g in order to find the flat region were the 
corrections are small enough, in order to obtain estimates of v with a smaller 
systematic error. Our final estimates are in reasonable agreement with the 
proposed value of u, the discrepancy being less than a few per cent. The 
worst cases are for A = 1 and A = 0.5. For instance for 8 = 1/3, A = 1 our 
data suggest v ~ 0.68,0.69 while the expected value is v = 2/3 ~ 0.667. 
Let us notice however that for A = 1 and A = 0.5 we expect the presence of 
logarithmic corrections, i.e. a behaviour of the form 

(R 2 e ) N = aN 2u \og^N(l + 0(l/logiV, log log N/ log N)) (48) 

The presence of logarithms makes the analysis very difficult. First of all 
in this case the convergence to the asymptotic regime is extremely slow. 
Moreover the presence of the term log^ 3 iV makes impossible an evaluation of 
v. Indeed as we use data with 200 < N < 8000, log^ N behaves, as far as the 
fit is concerned, approximately as jV~ a3/3 . Thus in a pure power-law fit one 
really measures 2z/ + 0.3/5. Since we don't have any theoretical knowledge of 
(3 it is thus impossible to draw any definite conclusion. 

To better understand the validity of our conjecture (|^) we have made 
two runs with higher statistics at A = 0.75 and A = 0.25 with 5 = 0. Each 
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Rl 


V 


Rl 


V 


N 


A = 0.25 


A = 


0.75 


500 


4862(10) 


0.7433(11) 


1874(4) 


0.6238(6) 


1000 


13554(31) 


0.7447(17) 


4449(10) 


0.6239(8) 


2000 


37927(91) 


0.7460(27) 


10556(25) 


0.6242(13) 


4000 


106704(277) 




25073(60) 




8000 


300005(874) 




59567(155) 




u(X) 




0.7500 




0.6250 



Table 2: Results for the runs with higher statistics and comparison with the 
expected value for the exponent v. In both cases 5 = and g = 1. 



data point corresponds here to 9 x 10 6 iterations and in order to avoid any 
initialization bias we have discarded the first 100 N iterations. The results 
are reported in Table |2[ A good agreement is seen although a systematic 
trend is visible in both cases. 

We have then checked if (Rl) obeys the scaling laws ( |30"D and fl3"9"l) . In 
Fig. p] we plot our estimates of (Rl) /N versus 

N l+&-\ for S < Q A = Q^QQ 

0.25, 0.50, 0.90 with A" > 200. The agreement seems quite good. However 
a more detailed examination of the scaling plot shows that the data points 
do not belong to a unique curve within error bars. Indeed one sees that 
different runs with the same values of g and 5 belong to distinct curves 
which approach each other only when N becomes large. This fact has to be 
expected. Indeed, as the analysis of the exponent v shows, for the values 
of A" which we are considering the corrections to scaling are still large and 
thus we expect analogously large violations to the scaling behaviour given by 
(j39|) . Analogously in Fig. ||we plot our estimates (Rl)/N 3 ^ 2 versus N 1 ^ 2+5 ^ x 
for S < and A = 0.75, 0.90 and A" > 200. Here again the agreement is only 
approximate. The situation is even worse for A = 1.00 and 5 > 0: in this 
case the points are scattered and no scaling can be observed. This can be 
explained by the presence of logarithmic terms, which break the scaling laws 
and make the approach to the asymptotic regime extremely slow. As a final 
check we have studied the behaviour of the universal ratio A, which, using 
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Figure 1: Scaling function for the end to end distance in the case 5 < 0. 
Small pentagons, rectangles, triangles and diamonds correspond respectively 
to the values A = 0.90, 0.50, 0.25 and 0. 
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Figure 2: Scaling function for the end to end distance in the case 5 > 0. 
Small pentagons and triangles correspond respectively to the values A = 0.90 
and 0.75 . 
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the scaling laws (|30| ) and (|3^) must have the form 

A = h x (gN 1+5 - x ) (49) 

for S < and 

A = hxigN 1 / 2 ^) (50) 

for 5 > 0. Moreover for x — > both functions must converge to the ORW- 
value 1/6. In Fig. |3| and f| we present the scaling plots for A in the two cases. 
The agreement is reasonable, although a closer inspection shows again the 
presence of systematic deviations. 



6 Conclusions 

In conclusion the results of our simulations are in good agreement with the 
theoretical arguments of Sect. 3, suggesting that our conjectured value for z/, 
if not exact, is certainly a very good approximation in the range of parameters 
we have examined. 

A Flory argument would predict for the exponent v 

» = 3 ~^^ (51) 
2 + D v ; 

It is interesting to remark that such a value is a good approximation, ac- 
cording to our analysis, only for A = 5 = 0. The deep reason for such an 
agreement escapes us. 
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Figure 3: Scaling function for the universal ratio A in the case 5 < 0. Small 
pentagons, rectangles, triangles and diamonds correspond respectively to the 
values A = 0.90, 0.50, 0.25 and 0. 
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Figure 4: Scaling function for the universal ratio A in the case 5 > 0. Small 
pentagons and triangles correspond respectively to the values A = 0.90 and 
0.75 . 
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g 


100 


200 


500 


1000 


2000 


4000 


8000 




A = 


-1 


10. 


207(1) 


414(2) 


1045(6) 


2071(12) 


4169(24) 


8284(48) 




-3/4 


2.0 


175(1) 


369(2) 


1016(5) 


2187(12) 


4702(26) 


10158(58) 




-3/4 


5.0 


249(1) 


541(3) 


1504(8) 


3268(18) 


7084(38) 


15342(86) 




-3/4 


8.0 


307(1) 


669(3) 


1874(10) 


4062(122) 


8876(50) 


19390(108) 




-1/2 


1.0 


207(1) 


477(3) 


1460(8) 


3475(19) 


8148(46) 


19273(108) 




-1/2 


3.0 


332(1) 


786(4) 


2480(13) 


5871(33) 


14000(80) 


33389(198) 




-1/2 


5.0 


417(1) 


1000(5) 


3164(18) 


7624(44) 


18227(110) 


42895(264) 




-1/4 


1.0 


341(1) 


876(5) 


3071(17) 


7980(46) 


20780(124) 


54207(350) 


140178(932) 


-1/4 


2.0 


462(2) 


1203(7) 


4285(25) 


11087(68) 


28817(188) 


75441(518) 




-1/4 


4.0 


612(2) 


1621(9) 


5815(36) 


15225(100) 


39867(274) 


104450(754) 






A = 0.25 


-3/4 


10. 


257(1) 


519(3) 


1318(8) 


2671(15) 


5372(31) 


10747(64) 




-1/2 


5. 


322(1) 


716(4) 


2088(12) 


4689(27) 


10379(61) 


23336(138) 




-1/2 


10. 


458(2) 


1052(6) 


3092(18) 


6964(41) 


15736(96) 


35485(214) 




-1/2 


34. 


726(3) 


1859(11) 


6034(39) 


14344(94) 


32933(222) 


74819(522) 




-3/8 


5. 


430(2) 


1035(6) 


3251(19) 


7746(47) 


18399(116) 


43618(270) 




-3/8 


13. 


658(3) 


1667(10) 


5513(35) 


13613(89) 


32933(222) 


79261(456) 




-1/4 


0.7 


212(1) 


507(3) 


1615(9) 


4001(23) 


9708(57) 


23969(140) 




-1/4 


4.7 


552(2) 


1405(8) 


4842(30) 


12376(81) 


31353(214) 


78563(550) 







1.0 


458(2) 


1259(7) 


4842(30) 


13613(89) 


37872(266) 


106280(806) 






A = 0.5 


-1/2 


22. 


529(2) 


1176(7) 


3174(20) 


6640(42) 


13593(49) 


27748(167) 




-1/4 


1.0 


193(1) 


428(2) 


1236(8) 


2778(17) 


6163(38) 


13877(86) 




-1/4 


3.3 


354(1) 


829(5) 


2525(16) 


5912(37) 


13593(49) 


31605(204) 


72165(472) 


-1/4 


6.0 


488(2) 


1168(7) 


3690(22) 


8784(58) 


20578(138) 


48194(330) 


112036(776) 





0.5 


238(1) 


592(3) 


2026(12) 


5198(32) 


13593(49) 


35624(232) 







1.0 


342(1) 


888(5) 


3174(20) 


8390(55) 


22322(151) 


59257(410) 







2.0 


502(2) 


1339(6) 


4925(23) 


13215(64) 


35925(190) 


96580(528) 


262693(1550) 





10. 


767(2) 


2133(12) 


8320(56) 


23177(170) 


64877(506) 


182090(1524) 


506839(4560) 




A = 0.9 


-1/10 


0.2 


115.4(0.6) 


237.1(1.6) 


608(4) 


1228(9) 


2528(20) 


5134(36) 


10465(74) 


-1/10 


0.5 


138.4(0.6) 


289.7(2.0) 


766(5) 


1594(11) 


3242(22) 


6743(46) 


13839(96) 


-1/10 


1.0 


174.1(0.8) 


373.2(2.5) 


1006(7) 


2119(14) 


4419(30) 


9287(65) 


19391(135) 



Table 3: Values of R 2 e for 5 < 0. All error bars are two standard deviations. 
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5 g 


100 


200 


500 


1000 


2000 


4000 


8000 












A = 


16 









1.0 


304 


;i) 


754(5) 


2582(16) 


6548(42) 


16710(113) 


42671(303) 







2.0 


448 


[2) 


1156(7) 


4080(27) 


10614(74) 


27816(205) 


71843(552) 







4.0 


620 


is) 


1653(10) 


6075(42) 


16173(118) 


43199(344) 


113130(980) 







5.0 


664 




1800(11) 


6672(48) 


17927(136) 


48107(390) 


129357(1110) 






A = 0.75 





0.5 


182 


;i) 


412(3) 


1236(8) 


2816(19) 


6571(44) 


15224(100) 







1.0 


253 


;i) 


599(4) 


1865(12) 


4473(30) 


10599(72) 


24984(174) 







5.0 


593 


is) 


1523(10) 


5253(39) 


13515(106) 


33763(288) 


84913(750) 




1/10 


0.5 


224 


;i) 


541(3) 


1784(12) 


4479(30) 


11058(76) 


28044(202) 




1/10 


1.0 


326 


;i) 


829(5) 


2845(19) 


7275(52) 


18615(178) 


47802(370) 




1/10 


2.0 


481 


[2) 


1262(8) 


4503(32) 


11705(90) 


30741(250) 


79394(680) 




1/4 


1.0 


479 


'2) 


1330(9) 


5117(38) 


14303(112) 


40132(342) 


111634(1042) 




1/4 


2.0 


642 


[2) 


1773(12) 


6842(50) 


19183(148) 


53193(450) 


147806(1320) 














A = 


19 









0.3 


136 


;i) 


290(2) 


789(5) 


1678(11) 


3605(25) 


7703(54) 


16455(114) 





0.5 


160 


;i) 


348(2) 


967(7) 


2127(15) 


4579(32) 


10022(68) 


21777 (145) 





1.0 


214 


;i) 


479(3) 


1400(9) 


3131(22) 


6943(51) 


15525(110) 


34318 (244) 


2/10 


0.2 


159 


;i) 


367(2) 


1141(8) 


2727(19) 


6593(46) 


16116(116) 


39180(288) 


2/10 


0.5 


240 


;i) 


596(4) 


1990(14) 


5003(37) 


12662(96) 


31804(252) 


80928 (675) 


2/10 


1.0 


353 


[2) 


911(6) 


3169(24) 


8116(64) 


21001(176) 


53719(482) 


135749 (1290) 












A = 


1 









0.1 


109.2 


;o.5) 


224.0(1.6) 


571(4) 


1167(8) 


2393(17) 


4861(34) 


9866(70) 





0.5 


148.5 


;o.7) 


318.8(2.2) 


857(6) 


1810(12) 


3824(27) 


8104(56) 


17021(122) 





1.0 


192.8 


;o.9) 


426.3(2.8) 


1190(8) 


2569(18) 


5541(39) 


11879(84) 


25536(186) 


1/3 


0.2 


187.7 


'0.9) 


462(3) 


1563(11) 


4016(30) 


10389(80) 


26944(220) 


70143(616) 


1/3 


0.5 


299.3 


:i-4) 


784(5) 


2829(22) 


7512(62) 


19567(172) 


51537(478) 


133100(1334) 


1/3 


1.00 


438 




1176(8) 


4245(34) 


11248(96) 


29304(277) 


77797(771) 


203511(2172) 


1/3 


2.15 


600 




1621(11) 


5906(46) 


15923(133) 


42261(385) 


112434(1100) 




1/2 


1.41 


686 


7) 


1907(18) 


7389(71) 


20805(210) 


55614(492) 


154556(1444) 


435115(4298) 


1/2 


1.79 


729 


;§) 


2038(25) 


7901(70) 


22182(208) 


58609(504) 


163330(1532) 


459835(4484) 


1/2 


2.83 


706 


;«) 


2000(23) 


7897(70) 


22483(208) 








2/3 


1.00 


655 




1820(12) 


7091(86) 


19756(185) 









Table 4: Values of B? e for 5 > 0. All error bars are two standard deviations. 
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$ g 


100 200 500 1000 2000 4000 


Vth 




A = 




— 1 10. 


5003(81 500014 2"! 4985(191 4999(301 4953(591 


5000 


-3/4 2.0 
—3/4 5.0 
-3/4 8.0 


0.5505(8) 0.5532(12) 0.5534(18) 0.5538(29) 0.5556(58) 
5588(81 5585(121 5584(181 5578(291 5575(561 
0.5617(8) 0.5616(12) 0.5621(18) 0.5638(28) 0.5637(57) 


0.5625 


— 1/2 1.0 
-1/2 3.0 
-1/2 5.0 


6141(81 6178(121 61QQ(181 617Q(281 6210(571 
6246(81 6255(121 6253(191 6269(301 6270(591 
0.6290(8) 0.6284(12) 0.6272(19) 0.6232(30) 0.6174(62) 


6250 


-1/4 1.0 
-1/4 2.0 
-1/4 4.0 


0.6869(7) 0.6884(10) 0.6895(14) 0.6894(20) 0.6886(32) 0.6854(67) 
0.6902(9) 0.6901(13) 0.6894(20) 0.6915(33) 0.6942(68) 
0.6969(9) 0.6952(14) 0.6945(22) 0.6946(35) 0.6948(72) 


0.6825 




A = 0.25 




-3/4 10. 


0.5068(8) 0.5062(12) 0.5046(19) 0.5022(30) 0.5002(60) 


0.5000 


-1/2 5.0 
-1/2 10. 
-1/2 34. 


0.5807(8) 0.5810(12) 0.5797(19) 0.5787(30) 0.5844(60) 
0.5897(8) 0.5872(13) 0.5869(19) 0.5873(31) 0.5866(62) 
6329(91 6177(141 6051(221 5958(351 5919(701 


0.5833 


-3/8 5.0 
-3/8 13. 


6263(81 6246(131 6243(201 6233(311 6226(641 
0.6513(9) 0.6453(14) 0.6408(21) 0.6354(35) 0.6335(70) 


6250 


-1/4 0.7 
-1/4 4.7 


0.6400(8) 0.6439(12) 0.6476(19) 0.6457(30) 0.6520(60) 
0.6732(9) 0.6723(14) 0.6703(21) 0.6667(34) 0.6626(71) 


0.6666 


1.0 


0.7378(9) 0.7404(14) 0.7423(22) 0.7410(36) 0.7443(75) 


0.7500 




A = 0.5 




-1/2 22. 


0.5382(8) 0.5277(12) 0.5212(20) 0.5158(32) 0.5147(53) 


0.5000 


-1/4 1.0 
-1/4 3.3 
-1/4 6.0 


0.5794(9) 0.5805(13) 0.5808(20) 0.5801(31) 0.5855(63) 

0.6076(7) 0.6061(10) 0.6047(15) 0.6026(20) 0.6030(26) 0.5956(66) 

0.6215(7) 0.6187(11) 0.6153(15) 0.6124(22) 0.6112(35) 0.6085(70) 


0.6250 


0.5 
1.0 
2.0 
10.0 


0.6774(8) 0.6835(12) 0.6893(19) 0.6942(33) 0.6950(54) 
0.6980(9) 0.7012(14) 0.7039(21) 0.7051(34) 0.7042(70) 
0.7139(6) 0.7153(8) 0.7171(12) 0.7183(17) 0.7174(29) 0.7218(59) 
0.7410(8) 0.7417(12) 0.7415(18) 0.7422(26) 0.7416(43) 0.7384(89) 


0.7500 




A = 0.9 




-1/10 0.2 
-1/10 0.5 
-1/10 1.0 


0.5142(8) 0.5135(12) 0.5138(16) 0.5149(23) 0.5124(38) 0.5137(72) 
0.5258(8) 0.5236(11) 0.5216(16) 0.5206(22) 0.5235(35) 0.5186(70) 
0.5381(8) 0.5352(11) 0.5334(16) 0.5326(22) 0.5334(35) 0.5311(72) 


0.5000 



Table 5: Values of z/ e // for 5 < 0, for various values of N cut . 
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5 g 


100 200 500 1000 2000 4000 


Vth 




A = 0.60 




1 
2 
4 
5 


0.6698(9) 0.6733(14) 0.6745(22) 0.6761(35) 0.6763(70) 
0.6886(9) 0.6898(15) 0.6903(23) 0.6899(37) 0.6845(77) 
0.7072(10) 0.7063(16) 0.7041(25) 0.7020(41) 0.6945(85) 
0.7147(10) 0.7134(16) 0.7128(25) 0.7128(41) 0.7135(85) 


0.7000 




A = 0.75 




0.5 
1.0 
5.0 


0.5996(9) 0.6023(14) 0.6045(21) 0.6086(34) 0.6061(68) 
0.6229(9) 0.6234(14) 0.6239(22) 0.6204(35) 0.6186(70) 
0.6745(10) 0.6718(16) 0.6687(26) 0.6628(43) 0.6653(89) 


0.6250 


1/10 0.5 
1/10 1.0 
1/10 2.0 


0.6537(9) 0.6582(14) 0.6612(22) 0.6614(35) 0.6713(72) 
0.6754(10) 0.6766(15) 0.6783(24) 0.6790(38) 0.6803(77) 
0.6927(10) 0.6918(16) 0.6908(25) 0.6907(41) 0.6844(85) 


0.6750 


1/4 1.0 
1/4 2.0 


0.7385(10) 0.7397(17) 0.7416(27) 0.7413(44) 0.7380(91) 
0.7375(10) 0.7386(16) 0.7387(26) 0.7364(42) 0.7372(89) 


0.7500 




A = 0.90 




0.3 
0.5 
1.0 


0.5470(8) 0.5476(11) 0.5481(15) 0.5487(22) 0.5476(35) 0.5475(70) 
0.5608(8) 0.5608(11) 0.5609(16) 0.5598(23) 0.5624(36) 0.5598(71) 
0.5798(8) 0.5788(11) 0.5770(16) 0.5762(23) 0.5763(36) 0.5722(72) 


0.5500 


2/10 0.2 
2/10 0.5 
2/10 1.0 


0.6271(8) 0.6336(11) 0.6382(16) 0.6411(23) 0.6428(37) 0.6408(74) 
0.6631(8) 0.6656(12) 0.6679(17) 0.6687(25) 0.6689(41) 0.6737(83) 
0.6803(9) 0.6793(13) 0.6787(19) 0.6778(28) 0.6733(46) 0.6687(94) 


0.6500 




A = 1 




0.1 
0.5 
1.0 


0.5141(8) 0.5136(12) 0.5139(16) 0.5131(23) 0.5109(36) 0.5106(72) 
0.5411(8) 0.5394(12) 0.5393(16) 0.5392(23) 0.5386(36) 0.5353(72) 
0.5579(8) 0.5546(11) 0.5529(16) 0.5520(23) 0.5511(37) 0.5521(73) 


0.5000 


1/3 0.20 
1/3 0.50 
1/3 1.00 
1/3 2.15 


0.6741(8) 0.6811(13) 0.6859(18) 0.6876(26) 0.6888(42) 0.6902(86) 
0.6972(9) 0.6969(14) 0.6949(20) 0.6921(29) 0.6918(48) 0.6844(99) 
0.7014(10) 0.6987(14) 0.6979(21) 0.6969(31) 0.6992(51) 0.6937(105) 
0.7099(11) 0.708(18) 0.7083(28) 0.7049(46) 0.7058(97) 


0.6666 


1/2 1.41 
1/2 1.79 
1/2 2.83 


0.7348(13) 0.7341(16) 0.7326(22) 0.7320(32) 0.7418(48) 0.7466(98) 
0.7332(13) 0.7315(17) 0.7302(21) 0.7300(31) 0.7428(47) 0.7467(98) 
0.7514(29) 0.7518(46) 0.7547(92) 


0.7500 


2/3 1.00 


0.7402(27) 0.7410(35) 0.7391(111) 


0.7500 



Table 6: Values of i/ e ff for 5 > 0, for various values of N cut . 
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